Search for the "ultimate state" in turbulent Rayleigh-Benard convection 
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Measurements of the Nusselt number Nu and of temperature variations AT;, in the bullc fluid are 
reported for turbulent Rayleigh-Benard convection of a cylindrical sample. They cover the Rayleigh- 
number range 10^ < iia < 3 x lO" using He (Prandtl number Pr = 0.67), N2 {Pr = 0.72) and SFe 
(Pr = 0.79 to 0.84) at pressures up to 15 bars and near-ambient temperatures. The sample had 
a height L — 2.24m and diameter D = 1.12m and was located in a new High-Pressure Convection 
Facility (HPCF) at the Max Planck Institute for Dynamics and Self-Organization in Gottingen, 
Germany. The data do not show the transition to an "ultimate regime" reported by Chavanne et 
al. and are consistent with the measurements of Niemela et al. 

PACS numbers: 47.27.te, 47.27.-i, 47.55. P- 



Turbulent convection in a fluid heated from below 
(Rayleigh-Benard convection or RBC) [H plays a major 
role in numerous natural processes. It occurs in Earth's 
outer core 0, atmosphere [4, 5], and oceans @, Q, 
and is found in the outer layer of the sun Q and in giant 
planets Q . The intensity of the driving by the thermal 
gradients usually is expressed by the dimensionless tem- 
perature difference known as the Rayleigh number Ra 
(to be defined explicitly below). Another important di- 
mensionless parameter is the ratio of viscous dissipation 
to thermal dissipation known as the Prandtl number Pr. 
The natural phenomena mentioned above generally in- 
volve Ra > 10^" and a wide range of Pr (see, for instance, 
Ref. [10]), whereas measurements in Earth-bound labo- 
ratories, with a few exceptions to be discussed below, had 
been limited to Ra < lO^^fH 

Of particular interest has been the global heat trans- 
port by the turbulent system, as expressed by the Nus- 
selt number Nu (the ratio of the effective conductivity 
Xeff of the convecting system to the conductivity A of 
the quiescent fluid). Extrapolations of Nu{Ra,Pr) to 
the geo- or astro-physically relevant ranges are in ques- 
tion because the basic physics involved in the turbu- 
lent flow is expected to change at some Pr-dependent 
Ra = Ra*{Pr). For Pr ~ 1 Ra* is estimated to be near 
3 X 10^^, and it is expected to increase with Pr approx- 
imately as Pr°'^.[lH Unfortunately, on this important 
issue two nominally equivalent sets of Nu measurements 
I3| 3 [I3| in the parameter ranges where Ra* might be 
found disagree with each other. One of them, by a group 
in Grenoble, France 12, 1^ (the "Grenoble" 




data) and 

shown in Fig. [T] as plusses, was interpreted as evidence 
for Ra* at an unexpectedly low value near 10^^ (this ev- 
idence will become more obvious in Fig. [3] below). The 
other, by a group in Oregon, 1^, (the "Oregon" data) 
shown in Fig. [T] as stars, did not show any transition all 
the way up to Ra ~ 10^*". Both experiments were done 
at approximately 5 K using helium near its critical point. 
Their relationships to Ra* are somewhat uncertain. Be- 



FIG. 1: The Nusselt number Nu as a function of the Rayleigh 
number Ra on logarithmic scales. Plusses (red online): from 
Chavanne et al. [T3|. Stars (purple online): from Niemela et 
al. [lj|] after a re-analysis reported in Ref. 16]. Solid dia- 
monds (red online): this work. He. Solid circles (blue online): 
this work, N2. Solid squares (black online): this work, SFe. 



cause of the proximity of the critical point Pr increased 
significantly as Ra increased beyond about 10^^ and thus 
the expected value of Ra* increased. 

Here we report new measurements, made at ambi- 
ent temperatures using compresses helium (He), nitro- 
gen (N2), and sulfur hexafluoride (SFg), over the range 
10^ < Ra ^ 3 X 10^*. They are shown as solid sym- 
bols in Fig. [TJ Our data fall into the narrow range 
0.67 < Pr ^ 0.84 characteristic of classical gases where 
theoretical estimates |ll| give Ra* ~ 2 x 10^'', albeit with 
a considerable uncertainty. Our measurements do not 
reveal the anticipated transition to an ultimate regime; 
they are roughly consistent with the Oregon data but are 
inconsistent with the Grenoble data. It remains unclear 
to us what caused the transition revealed by the Grenoble 
results. 

For modest, laboratory-accessible, Ra the heat trans- 
port is controlled essentially by thin thermal boundary 
layers (BLs) just above (below) the bottom (top) plate 
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of the convection cell. For Ra > Ra* it is expected 
that a large-scale circulation (LSC) in the sample in- 
terior will apply sufficient shear to the BLs to cause 
them to undergo a turbulent transition. This transition 
was predicted to occur when the shear Reynolds num- 
ber Res, based on the BL thickness, exceeds a critical 
value i?e* ~ 420. [ll| It is expected to lead to a new tur- 
bulent convective state in the sample interior that was 
analyzed theoretically by Kraichnan [l3| and that is of- 
ten referred to as the "ultimate state" because it should 
prevail asymptotically up to infinite Ra. Clearly the ex- 
perimental discovery and study of the ultimate state is 
one of the major challenges to experiment in this field, 
both because of its fundamental interest and because of 
its geo/astrophysical relevance. This is the reason why 
understanding the difference between the Oregon and the 
Grenoble data is a major issue in the field of turbulent 
convection. 

Before proceeding to a discussion of our results, it is 
useful to re-examine the prior large- i?a measurements in 
more detail. Over 30 years ago it was recognized that the 
study of RBC using low-temperature helium has unique 
advantages because of the unusual properties of this fluid. 
[liL ligl . liol l This was exploited in early experiments by 
a group in Chicago who reached Ra ~ 6 x 10^^ [^ij us- 
ing a cylindrical sample of aspect ratio F = D/L = 0.5 
{D is the diameter and L the height) by approaching 
the critical region of Helium where Ra tends to be- 
come large; but that work did not reveal any transition 
in the heat transport. Later the Grenoble group 
carried out similar measurements, also using a cylinder 
with F = 0.5, and tentatively identified a transition at 
Ra* ~ 10^^. Below that point Nu roughly followed an ef- 
fective power law Nu = Ra'^''^^ , with 7e// ~ 0.3. Above 
Ra* Jeff approached values near 0.4 and was heading 
toward the asymptotic value 0.5 predicted by Kraichnan. 
Soon thereafter the same group published a Physical Re- 
view Letter entitled "Observation of the ultimate regime 
in Rayleigh-Benard convection" . [13] A number of other 
papers by this group followed with more detail and data 
up to Ra ~ 10^^, and we show the results reported in 
Ref. [15] in Fig. [T] as plusses. It remained unexplained 
why the Chicago group did not find the transition seen 
by the Grenoble group. 

A few years later the Oregon group extended the he- 
lium measurements to theretofore unprecedented values 
of Ra as large as 10^^ by building a cylindrical convection 
cell with a height L = 1.0 m and a diameter D = 0.5 m 
and operating near the critical point of helium. 43] Their 
data[l6| are shown as stars in Fig.[T]and do not reveal any 
transition. They can be described over the wide range of 
Ra by a power law with Jeff — 0.32 with no significant 
tendency for Jeff to increase with Ra. 

We note that the Chicago, Grenoble, and Oregon 
experiments were conducted at cryogenic temperatures 
where experimental difficulties are quite severe and where 



on average Pr tends to increased with Ra. It was very de- 
sirable to carry out an experiment at near-ambient tem- 
peratures using more conventional, and above all differ- 
ent, experimental techniques and fluids with essentially 
constant Pr that could also achieve the large- i?a values 
deemed necessary to reach Ra* . To satisfy those criteria, 
we used a very large pressure vessel at the Max Planck In- 
stitute for Dynamics and Self-Organization in Gottingen, 
Germany. It is a cylinder of diameter 2.5 m and length 
5.5 m, with its axis horizontal, and with a turret above it 
that extends the height to 4 m over a diameter of 1.5 m. 
Because of its suggestive shape, this vessel has become 
known as the "Uboot of Gottingen" . It has an approxi- 
mate volume of 25 m'^ and can be used up to pressures of 
15 bars. In the section containing the turret we placed a 
RBC sample-cell with L = 2.24 m and = 1.12 m (the 
"High Pressure Convection Facility" or HPCF), yielding 
F = 0.500. It had top and bottom plates made of Alu- 
minum and a Plexiglas side wall of thickness 9 mm. The 
top plate extended from the top a short distance into the 
side wall; it was water cooled with flow through quadru- 
ple spiral channels. The bottom plate was a composite 
consisting of two aluminum plates with a Plexiglas plate 
of thickness 5 mm between them, all glued together with 
Stycast 1266 epoxy. The composite extended from the 
bottom into the side wall. Joule heating was applied 
with a heater that was uniformly distributed over the 
bottom surface of the bottom plate. Aside from that 
total heat current Q, the heat current Qs entering the 
sample could be inferred from the temperature difference 
across the composite. 23|, |2J] Under most conditions Qs 
was equal to Q, indicating no significant parasitic heat 
losses from the bottom-plate heater. In all cases Qs was 
used to compute Nu. Various thermal shields, regulated 
at appropriate temperatures to prevent heat losses, had 
been installed as described for instance in Ref. [1^. All 
spaces outside of the sample and up to a diameter of 1.4 m 
were occupied by low-density open-pore foam. Estimates 
indicate that side-wall heat-losses were negligible. Each 
data point was derived from a time series of temperature 
readings at time intervals of about 5 sec and spanning 
about one day. The first half of each series was discarded 
to avoid transients and the remainder was time-averaged 
before computing Ra and Nu. The data were corrected 
for the finite plate conductivity 25, [2^; this correction 
generally was less than a few percent. 

The Rayleigh number is given by Ra = PgL^AT/ (kz/), 
and the Prandtl number is Pr — v j k. Here the thermal 
expansion coefficient /?, the thermal diffusivity k, and the 
kinematic viscosity v were evaluated at the mean temper- 
ature r„i = (Tj + Th)/2 (Tt and Tf, are the temperatures 
at the top and bottom of the fluid respectively). The 
gravitational acceleration is and AT = Ti, — Tf. The 
Nusselt number is Nu = QsL/{AATX) where A is the 
sample cross-sectional area and A is the thermal conduc- 
tivity. The fluid properties were compiled previously 27 ] 



3 



10 



*» 



+ ^ 



+ + + + -H- 



i» -t 1 i i »bjiwtwN i i 



109 



10^ 



10 



13 



10 



15 



Ra 



FIG. 2: The Prandtl number as a function of the Rayleigh 
number corresponding to the various data sets. The symbols 
are as in Fig. [1] The dotted hne is the estimate of Pr{Ra*) 
from Ref. [ull for F = 1. 



from numerous papers in the literature. 

Measurements of Nu{Ra) were made with He {Pr — 
0.67), N2 {Pr = 0.72), and two separate gas fiUings (run 
1 and run 2) of SFg (Pr = 0.79 to 84). The data span 
the ranges 10^ < Ra < 3 x 10^"'. The results are shown in 
Fig. [1] as solid symbols. Figure [2] shows the points in the 
Ra — Pr parameter space where the various data were 
taken. Even though at large Ra the Prandtl number is 
much lower (and thus Ra* is much lower) than it is for 
the Oregon or Grenoble data, our data in Fig. [T] reveal 
no transition to an ultimate regime. The dotted line in 
Fig. [2] is an estimate of Ra* (Pr) for cylindrical samples 
with r = 1 which may serve as a rough guide of where to 
expect the transition also for F — 0.5. We see that the 
Oregon data and our new results exceed that estimate 
by about the same amount even though the Oregon data 
extend to much larger Ra. The Grenoble data, on the 
other hand, do not come near the estimate; thus it seems 
unlikely that the transition near Ra = 10" which they 
reveal is associated with a shear-induced turbulence tran- 
sition in the BLs. 

In Fig. [3] we show the data for Nu in more detail by 
dividing out the approximate Ra dependence i?a°'^. For 
Ra ^ 10^^ there is general consistency of the dependence 
of Nu on Ra between all three data sets. Our He data are 
about 6% below the Oregon and Grenoble data and about 
4% below our own N2 data. Possible causes for these 
differences may be systematic errors in the calibration 
of the thermal resistance of the bottom-plate composite, 
systematic errors in the thermodynamic and transport 
properties used to compute Ra and Nu for the various 
fluids, as well as an expected increase of Nu with Pr at 
constant Ra for Ra ^ 1 . For this "reduced" Nusselt num- 
ber the break in the Grenoble data at Ra ~ 10^^ is more 
apparent than it was in Fig.[T] Clearly this transition to a 
significantly larger effective exponent, near 7e// = 0.39, 
is not present in the Oregon data; nor can it be found 
in the Gottingen data. In general our data are slightly 
lower that the Oregon data, and this difference increases 
as Ra increases. In part this may again be due to sys- 
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FIG. 3: The reduced Nusselt number Nu/Ra^'^ as a function 
of the Raylei gh n umber Ra. Plusses (red online) : from Cha- 
vanne et al. [l5| . Stars purple online): from Niemela et al. 
[T^ after a re-analysis reported in Ref. [l^. Solid diamonds 
(red online): this work. He, P = 4.3 bars. Solid squares 
(black online): SFe, run 1. Solid circles (red online): SFs, 
run 2. Open symbols: This work, N2. Circles (purple online): 
P — 2 bars. Squares (red online): P = 6 bars. Down pointing 
triangles (blue online): P = 10 bars. Up pointing triangles 
(purple online): P = 15 bars. The vertical dotted line shows 
the location of a transition, see text. The area in the dashed 
rectangle is shown enlarged in Fig. 3^. 
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FIG. 4: Details from SFg, run 1, of (a) the reduced Nusselt 
number Nu/Ra'^''^ and (b) the temperature drop 2{{Tt) — 
{Tt})/AT in the bulk near the side wall. Open (solid) circles: 
P ~ 8 bars, decreasing (increasing) pressure. Solid squares: 
15 ^ P ^ 9 bars. The vertical dotted lines show the location 
of a transition, see text. 



tematic errors in the fluid properties that were used in 
the analyses, especially near the critical point of helium; 
but at the larger Ra a dependence of the Oregon Nu 
data on Pr may also be contributing. As illustrated in 
Fig. [21 Pr starts to increase for the Oregon data as Ra 
exceeds about 10^^, whereas for our data Pr is nearly 
independent of Ra. 

A notable feature of the Goettingen data is a sud- 
den change of the dependence of Nu upon Ra near 
Ra = 6.7 X 10^^. Hardly noticeable in Fig. [I] it becomes 
apparent with the much higher resolution of Fig. [3] where 
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it is indicated by a vertical dotted line. One data point, 
per chance taken precisely at that transition, yielded a 
Nu value 12% higher than the neighboring ones (see the 
isolated solid square in Fig. [3]). To look at this feature 
more closely, we show in Fig.|3^ the data from run 1, SFe, 
in the dashed rectangle of Fig. [3] on an enlarged scale. 
There the circles are for P ~ 8 bars and the squares are 
for higher pressures. The open (solid) circles were taken 
with decreasing (increasing) Ra (and thus AT). Again 
the transition is indicated by the vertical dotted line. We 
believe that it reflects a change in the structure of the 
LSC, perhaps from a single roll to a two-roll structure as 
found earlier in direct numerical simulations 2^ 2^ and 
experiments [30| . but more work is required to be more 
specific. 

An interesting and instructive feature of this system 
is the vertical temperature drop AT/, across the bulk of 
the fluid. [3l| Whereas the temperature drop ATbl across 
the two boundary layers dominates AT = ATbl + ATf,, 
values of order 0.1 have been found for AT^/AT in exper- 
iments with water, Ra near 10^, and near the side wall 
where plumes travel from one plate to the other (any 
temperature drop along the sample axis where plumes 
are scarce was found to be considerably smaller). In the 
range Ra < Ra* the gradient in the bulk is expected 
to remain small; but as Ra* is exceeded, the BLs are 
expected to disappear and thus most of the tempera- 
ture drop should occur across the bulk. Thus we an- 
ticipate a dramatic change of AT;, as Ra* is exceeded. 
With this in mind, the azimuthal and temporal mean 
temperatures (T^) and (T() along the side wall at the 
vertical positions — L/4 and T/4 respectively were deter- 
mined with thermometers imbedded in the side walls as 
described for instance in Ref. [sij. Those data yield the 
estimates AT^ ~ 2((Th) - (Tt)) shown in Fig. [Id. There 
is some potentially interesting structure in the Ra depen- 
dence of this quantity, especially near the transition at 
Ra — 6.8 X 10^^ which is indicated by the vertical dotted 
line; but the overall value of AT^/AT remains modest 
at all Ra and is comparable to similar results reported 
elsewhere. fsij The results imply that ATbl » AT^, 
consistent with Ra < Ra*, for our entire parameter 
range. 

In this Letter we report new measurements of the 
Nusselt number, and of the bulk temperature gradient 
near the side wall, for turbulent Rayleigh-Benard con- 
vection in a cylindrical sample of aspect ratio 0.5 over 
the Rayleigh-number range 10^ ^ Ra ^ 3 x 10^^ and for 
Prandtl numbers close to 0.8. The results do not reveal 
any evidence for a transition to the "ultimate" regime 
expected in the large- i?a limit. 
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